Critical speed-up vs critical slow-down: a new kind of relaxation oscillation with application to 

stick-slip phenomena 
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The equations for the sliding of a single block driven by an elastic force show numerically a fast and a slow 
step in their dynamics when a dimensionless parameter is very large, a limit pertinent for many applications. An 
asymptotic analysis of the solutions explains well the two sharply different steps of the stick-slip dynamics. The 
stick (slow) part takes place along a slow manifold in the phase space. But, in contrast with standard relaxation 
dynamics (of van der Pol type), the slow manifold is always formally attracting and the transition from slow 
to fast dynamics occurs because the slow dynamics has a finite time singularity breaking the assumption of 
slowness. This makes a new kind of relaxation oscillation. We show that the response of the stick-slip system to 
an external noise displays a progressive speed-up before the transition, in contrast with the well known critical 
slowing-down observed in the standard case. 



INTRODUCTION 

Many physical systems show relaxation oscillations (RO) 
with phases of slow motion interrupted by fast drifts. Exam- 
ples of RO in the real world are found in electro-mechanical 
devices yj], in purely electrical devices yfl and in self- 
oscillating circuits ||4[]. Besides these saddle-folding 131] sys- 
tems, it has been suggested by Brace and Byerlee in 1938 151] 
that earthquakes are RO of the stick-slip kind, appearing when 
two solids slide on each other Here we compare the mecha- 
nism of RO in stick-slip models and in saddle-folding ones for 
which a prototype dealing with plane flow RO is the slightly 
extended van der Pol equation in the limit /3 large, written as 
follows in its Lienard form 



y = X + a 



(1) 



In this limit, the slow dynamics is a drift along the slow 
manifold (SM) of Cartesian equation y ~ x — ^ interrupted 
at places (the folding spots) where this manifold is unstable 
against transverse perturbations. There are then two widely 
different time scales, the short one for the jumps away from 
the SM and the long one for the slow drift. In this class of 
models, the transition from fast to slow motion occurs by dy- 
namical saddle-node bifurcation ]jg]. Near the saddle-node 
transition (ie. at .t = ±1 and y = ±2/3) the loss of stability 
of the SM is generically the same as the one of a particle mov- 
ing with friction in a potential changing slowly as a function 
of time, a model recently investigated in view of predicting 
catastrophes l^. 

This picture of RO in the strongly nonlinear limit with a 
slow manifold and quick jumps outside of this manifold seems 
to make the currently accepted "paradigm" for RO in dynam- 
ical systems. However, by studying the equations pertinent 
for the stick-slip dynamics we have found another scenario of 
RO, which behaves differently near the transition, with possi- 
bly important consequences on precursor phenomena. 

Stick-slip phenomena occurring in solid friction are ubiq- 
uitous in real world, for example in geophysics and engineer- 



ing. Earthquake rupture have been attributed to stick-slip in- 
stability observed in laboratory experiments [5^, [^ and ][^. It 
should be a good example of RO because the two phases (stick 
and slip ) take place with widely different speed. Mathemat- 
ical models for solid friction, as proposed by Dieterich llOll 
and by Ruina ]|8|], make a fair representation of this physics, 
their numerical solutions showing RO with two widely sepa- 
rated time scales. In the present paper we show that in this 
family of models the transition from slow to fast does not 
occur because the SM becomes repulsive, contrary to what 
happens in the strongly nonlinear regime of van der Pol-type 
models. Here the SM is attractive in the whole space. We 
show that the transition occurs because the trajectory on this 
manifold goes to infinity in finite time. Before the critical 
time tc the adiabatic approximation becomes invalid, and the 
real trajectory escapes from the SM. Therefore the stick-slip 
model of Dieterich-Ruina-Rice (DRR) makes a new class of 
RO where the slow trajectory begins to accelerate a long time 
before tc, whereas it slows down shortly before tc in the case 
of equation ([T]i. It follows that the statistical properties of the 
responses to an external noise are completely different for the 
two classes of models. For stick-slip models the correlation 
time of the response decreases before tc (speed-up), whereas 
it increases (slowing-down) in the case of the standard saddle- 
node transition ]|7|] . 



DRR MODEL EQUATIONS 

In the DRR models a "state" variable, denoted as 9 below, 
describes the physical state of the surfaces sliding on each 
other (rugosities, asperities, etc.) and how it influences the 
friction. This state variable has its own dynamics, in agree- 
ment with the observations of time dependent solid/solid fric- 
tion. Those models have been used to analyze the sliding of a 
block on a flat solid, the block of mass M being driven by an 
external spring with the other end moving at constant speed 
vq, this being possibly analogous to fault slip on the bound- 
aries of Earth's crustal plate responsible of earthquakes 1 1 ill . 
We shall consider below the single-block problem. 



We introduce first the DRR set of equations of motion of a 
single sliding block coupled with Dieterich- Ruina rate-and- 
state dependent friction 181], 1 1111 and discuss their solution in 
the limit where the inertia of the block is very small, the limit 
of the stick-slip sliding. This set of ODE's (ordinary differen- 
tial equations) read 



and 



V - Wo, 



V = -— {ku + 6 + Ahv{v/vi)) , 



and 



D, 



{e + B\TL{v/vi)). 



(2) 



(3) 



(4) 



The first equation relates u, the position of the block, to its 
speed V. The constant uo is the difference of speed between 
the point holding the spring and the surface the block is slid- 
ing on. The force due to the external spring is proportional to 
u, with a spring constant k, and M is the mass of the block. 
Moreover 9, a contribution to the friction force, depends on 
the state of the surfaces facing each other, a function of the 
history of sliding, as given by the solution of the equation 
^. The intermediate equation (|3]l is a way of writing the 
dynamics of the block under the effect of the pulling force 
proportional to u, of the friction, plus the state dependent part 
proportional to 6. The quantity {0 + A\n{v/vi)) on its right- 
hand side is the friction force plus a constant, absorbed into 
a constant displacement added to u. Even after addition of a 
constant friction (independent on v), the DRR equations can- 
not be valid for all values of (m, 0, v). Because A is positive, 
the friction force is dominated at low speeds v by the loga- 
rithm, and so becomes very large positive, in the direction of 
V, a negative friction which is clearly nonphysical. Because 
the added constant is independent on the parameters used to 
write the equations above, one cannot tell, without knowing 
what is its value in a particular application, if the non physi- 
cal regime of negative friction is reached or not for this case. 
We assume that, in the regime we consider, negative friction 
never happens. Note that in case of very low speed, Lapusta et 
al. lll2ll propose to change the above relation v/vi — cxp(r/) 
between the velocity v and the friction tj, by the relation 
v/vi = 2 sinh(T/), that changes ln(?;) into \n{v + yjv^ + 1). 
The velocity v\ and the coefficient A, B and Dc are phe- 
nomenological quantities derived in principle from experi- 
ments. One can absorb the ratio vq/vi by adding another 
constant to u. Using the scaled quantities v/vq, 9/A,u/Dc, 
vot/Dc, one transforms the above equations into the mathe- 
matically convenient form lll3ll 



u — V — 1, 



= --f^[u+-{e + \n{v)) 



(5) 



(6) 



'+(! + £) lll(f )) 



(7) 



This version of the DRR equations keeps three dimension- 
less parameters, 7, e and ^, related to the quantities M, vq 
and to the other phenomenological parameters A, B and Dc 
derived in principle from experiments 181], Ml If . The geophys- 
ical literature gives e and ^ of order one, and 7 large. To give 
an order of magnitude, 7 is typically of order 10 for ice flows 
in Antarctica ill4ll . where the sliding phase lasts about 20 min- 
utes once a day, while for major earthquakes 7 ^ 10^ ^ with 
sliding phases of typically 5 seconds occurring once every 200 
years. 

When the parameter 7 becomes very large and I- is above 
unity, the fix point (u = 6' = 0, u = 1) of equations (|5]l- 
(|7]i undergoes a Hopf bifurcation, leading to a stable periodic 
solution, see figure ^a). At threshold the period is 



7!' 



(8) 



The period t increases noticeably with the control parameter, 
as illustrated by curves (b-c) in Figure ([T]) drawn for ^ ~ 0.7, 
it reaches a linear dependence with respect to e for given 7 , 
and increases nearly as 5 log(7) for given e. 

In the case of earthquakes, the typical parameter values ]|9t] 
Vq = 30 mm/year, Dc/vq — 2.7 years, ^ ~ 0.8, e = 1,7 = 
10* lead to a period equal to 75 years, or r = 28 in units 
of equations ©-(JTll. More generally, considering earthquakes 
separated by 30 to 200 years, it gives the dimensionless period 



10-70, 



(9) 



which is the range of period values considered below in our 
numerical calculations. 

u, V y,e 




Log7 



FIG. 1: (a) Periodic solution of equations iSTfl for ^ — 10^,^ = 0.7, 
e = 0.72. (b-c) period of the limit cycle for ^ — 0.7, in (b) as 
function of e for 7 = 10' ; in (c) as function of logjo (7) for e — 0.72. 

The striking point is the stick-slip character of the solutions 
with almost motionless long intervals (stick step) interrupted 
by short bursts of rapid displacement (slip step). In this limit 



the stiffness of the DRR equations makes them hard to solve 
numerically, so that most numerical simulations are done with 
7 of order unity, unrealistic for earthquake modeling. There- 
fore is is important to carry an asymptotic analysis of the so- 
lutions of this set in the large 7 limit, what is done below. 

Moreover we check our theoretical results by a systematic 
numerical investigation of the solutions for 7 in the range 
10-10"*, plus some calculations for larger values of 7 ( un- 
til lO'^). This was performed on Mathematica, by the "Stiff- 
nessSwitching" method using a pair of extrapolation methods 
as the default. The nonstiff solver uses the "ExplicitModified- 
Midpoint" base method, the stiff solver uses the "Linearlylm- 
plicitEuler" base method. 



SLOW AND FAST REGIMES 

stick step 

Consider first the stick regime, where the acceleration of 
the motion is very small. 



w<7 



(10) 



although there is always some sliding as illustrated in figure 
dlJa). The time duration of this step is very close to the period 
T of the limit cycle. 

The velocity v and the acceleration v are smaller than one 
in the main part of the sticking stage, they both becomes unity 
at time nearly equal to tc — 1. This defines the intermediate 
time scale to, the time interval separating the instant where 
v{t) = 1 and the catastrophe time tc where the velocity is 
maximum. For the whole range of parameters that we investi- 
gated numerically, we have found the relation 



to = l, 



(11) 



within the units of equations (|5]l-(|2]l. After reaching the value 
V = l,and V = 1 the velocity and the acceleration increase. At 
the end of the stick step, in the range defined by 1 ^ w ^ 7 
and equation (fTOl i. the solution is analytically tractable, see 
the section "matching slow and fast solutions" below. Note 
that in the case of earthquake, the intermediate time scale 
corresponds to few years (2.7 years in [91), that makes the 
DRR model completely different from the saddle-node model 
where the intermediate stage lasts a few hours only |7[1 . 

During the whole sticking episode defined by equation ( fTOl l. 
the coefficient of 7^ on the right-hand side of equation (|6j can 
be set to zero, that gives the relation 



^u + 9 + \n{v) = 0, 



(12) 



= e- («"+'') in the 



defining a surface S of equation v{u 
phase space, (u, v, 9). Putting the relation (fTZt in equations 
(|5j and (|7J, one get a set of two first order ODE's for u and 6 
in which the large parameter 7 has disappeared. 



ii = e-(«"+e) - 1 

^ = e-(«"+«)[e0 + (l + e)^u] 



(13) 



defining a 2D flow (two dimensional). 

From the point of view of the 3D dynamical system we 
started from, the 2D flow defined by equations ( fT3] l takes place 
on the surface (manifold) 5, a reduction of the original 3D flow 
done by Gu et al. jlSil in the same limit. The "stick" phase of 
the 3D dynamics fulfills equations (fTOl i. and (fT2]i-(fT3Tl. 

Let us show that in the 3D flow, the SM is attracting. The 
surface S splits the 3D phase space {u,v,d) into two sub- 
spaces, the one for which (^it + 9 + \n{v)) is positive and the 
one for which it is negative. In the large 7 limit, the dominant 
component of the time derivative of the 3D vector (w, v, 9) is 
V. The negative sign on the r.h.s. of equation (|6]l, ensures the 
that the SM is attractive from both subspaces, because taking 
V very large positive one sees that v, the dominant component 
of the "velocity", defined by the set Q, (HI and ^ is directed 
toward v negative, and conversely for {^u + 9 + ln(w)) nega- 
tive. This shows that the slow manifold is everywhere attract- 
ing in this limit 7 large. Therefore the classical "scenario" 
for van der Pol-type RO with an ejection out of a repelling 
SM cannot hold for DRR equations. We show later that in the 
DRR case, the 3D flow escapes from S because the 2D flow 
trajectory goes to infinity m finite time, see Figures (O. 



Slip step 

During the slip, both v and v are "large", of order 7 and 
7^ respectively. Neglecting (—1) on the right-hand side of 
equation (|5]l allows to eliminate the large parameter 7 from the 
equations by rescaling and addition. This is done by putting 
V = -fV(T), t = r/7, u = U{T) + Cu ln(7) and 9^Q{T) + 
cg In (7). The quantities denoted as T, V, 8 and U are of order 
1 with respect to 7 as well as the constants cq^Cu- Setting to 
zero the part independent of T in the equations derived from 
^ and (|7| one finds two algebraic equations for eg and c„ 
with the solution 






-(1 + e) 



(14) 



?' 



The functions U{T),V{T) and 8(T) are solutions of a set of 
three ODE's free of the large parameter 7 



Vt -- 



where Xt is for ^. 



V 



U+^{e + ln{V)) 
V^ [6 + (1 + e) \n{V) 



(15) 



This predicts that the slip lasts a (short) time St^qk of order 
1/7 although the sticking lasts a time independent on 7. This 
scaling law is in excellent agreement with the numerical study 
of the full DRR equations for 7 large. We observed that, as a 
function of t, v{t) peaks with a time half- width 

Stcqk - 2.5/7, (16) 

in a wide range of parameters, as illustrated in Figure (O. 
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FIG. 2: Fast step. Part of the solution shown in Figure ([TJ, close to 
the catastrophe, in units of equations l ll5b 



MATCHING SLOW AND FAST SOLUTIONS 

There remains to interpolate between the stick and shp parts 
of the dynamics. To the best of our knowledge, this matching 
has not been done before for DRR-like equations, in relation 
with the idea of a finite time singularity of the slow dynamics. 
According to matching theory, the solution in the two interpo- 
lation domains (before and after the catastrophe) has to satisfy 
the set of equations common to the two phases, slow and fast. 



U = V, 



(17) 



together with equations (|7]l and ( fT2] i. The integral curves of 
dTji and ([TtT i are solutions of the single linear ODE 



that gives 



du 



e{u) = 



= e9 + {l + e)^u, 



(1 + e)? 



(1 



(18) 



(19) 



with c arbitrary constant. In our problem this expression 
agrees well with the numerics, with c calculated from a partic- 
ular value of (u, ff) belonging to the matching region, see the 
discussion below concerning the figures (|3}a-b) and (|2|i. The 
general time dependent problem in the matching region may 
be written as m = exp($(w)) with 



$(ii) 



—u 
e 



(1 + ^K 



(20) 



The solution t ~ to — J,^ , du'Gxp(— $(w')) yields a 
singularity at finite time in the forward time direction if 
/ du' exp(— <&(u')) is a converging integral, which requires 
c to be negative. This singularity occurs at the end of the 
stick interval, as illustrated in figure (|3}a). However we must 
note that the solution ( |20l i was derived within the adiabatic 
approximation (fTOl i, which limits its validity to a domain we 
are going to precise. On the surface S, the time derivative of 
V = exp{-^u - 9) writes i = e^^lC^+e) [^ - e0 - {1 + e)v] 
by using equationdTTb and the second equation ( fT3] l. There- 
fore in the space (u, 9) the straight line of Cartesian equation 




(a) ll.b^O n.635 17,640 17.645 17650 17.655 17.660 




FIG. 3; (a) v{t)/^ for the 2D and 3D flows (blue and purple curves 
respectively) (b) Phase space (u, 9). The 3D limit cycle (closed 
dashed blue curve) is superimposed with the 2D flow (solid red por- 
tion) on the SM. The orange and blue portions correspond to the two 
branches of matching solutions, equation l ll9t . 



splits the SM into a part such that trajectories crossing this 
line tends to still lower values of ii, and so follow even more 
closely the SM, although when the crossing is in the other di- 
rection, i) tends to increase, as well as v itself (always positive) 
and so ends up with the finite time divergence just described. 
A rough estimate of condition ( fTOl ) limits the 2D flow validity 
to the domain located above the line 



(u) = -Cu - ln(7). 



(22) 



e9{u) ^ {1 + €)u - ^. 



(21) 



These results are illustrated in figure ([3}b). In this figure 
the closed curve (blue-dashed) displays the limit cycle (3D 
flow) evolving counterclockwise in the phase space {u,9). Su- 
perimposed onto this curve, we plot the stick episode (solid 
red portion), and the matching analytical solutions iT% (or- 
ange and green curves). The solid straight line corresponds to 
the zero acceleration relation ( 1211 1. joining the two extrema of 
the closed curve 9{u). The adiabatic approximation domain 
lies above the dashed straight line, equation ( |22] |. The two 
portions of analytical solutions fit well the two intermediate 
regimes, they are calculated with two distinct values of c (de- 
rived from equation ( fT9b and using values of {u, 9) before and 
after the catastrophe respectively). The orange curve which 
belongs to the SM, matches the stick regime and the fast one, 
then it diverges on the surface 5*, while the limit cycle escapes 
from S. This happens when the adiabatic approximation fails. 
The escape from the SM of the 3D flow in the fast regime is 
visible on figure (|4]i. 

The physics of the oscillation can be understood as follows. 
The acceleration of the block at the end of the stick phase 
is unstable because it lowers the friction and so feeds itself. 




V/7 



FIG. 4: 3D view of the trajectory (solid red curve) leaving and land- 
ing out and in the SM (yellow-green surface), the black part diverging 
on S displays the 2D flow in its "non slow" regime. 



The transition to the "sliding" fast regime occurs when the 
acceleration becomes big enough to make the inertia of the 
block relevant. In this fast regime the dynamical system gets 
one more dimension (in other terms the effect of inertia in- 
creases by one the number of coupled ODE's). If one neglects 
in this fast regime the friction, the dynamics is the one of a har- 
monic oscillator making a half swing while its velocity keeps 
the same sign. When this oscillator has slowed down enough 
to yield back its energy to the potential energy of the spring, 
its velocity returns to small values and friction becomes so 
large that the SM is reached again. This could explain that, 
in earthquakes, a finite fraction of the initial elastic energy is 
not dissipated in the event itself, because part of it remains as 
potential energy after the large scale pendulum-like motion. 
In this model, the spring is under dilative strain before the slip 
episode and under compressive strain at the end of it. 



LYAPUNOV ANALYSIS ALONG THE TRAJECTORY 

We have noted above that the SM is everywhere attractive, 
taking argument of the negative sign in the r.h.s. of equation 
(|6]l. This statement is invalid for the close vicinity of the SM, 
more precisely at distance of order 1/7^ from the SM. In the 
very slow part of the trajectory, we may investigate how attrac- 
tive is the SM by performing the linear stability analysis of the 
flow. Along a peculiar trajectory {uo{t) ,hi{vo{t)) , 9o{t)) the 
jacobian matrix is 
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7- 


i'vo 


? 


A{9o,vo) 


-Vo 



(23) 



A^ + (^ + vo)X^ + 7^^T^ + 1)A + l^vo = 0, (25) 
In the large 7 limit, the two eignevalues are 



A(i) = -Bit) ± ^B{ty-^vo{ty 



with 



B{t) = -voil-Aieo,vo) 



-0- 



(26) 



(27) 



The three Lyapunov exponents are drawn as functions of time 
in figure (|5]l, for the flow coiTesponding to figure ([T}a), or ([3} 
b) or (HJi. One of the eigenvalue (curve a) is negative all along 
the trajectory. The two other eigenvalues are real in the first 
part of the trajectory, and become complex conjugate at time 
about tc — 4.6, see figure (b). The real parts are the red and 
blue curves in (b) , the imaginary parts are drawn in figure 
(c). The real part crosses zero at time t ^ tc — 2.3, with the 
eigenfrequency lo of order unity. This occurs inside the "very 
slow" stick regime ending near t ^ tc — to where v = 1. At 
this time, the Lyapunov analysis becomes invalid, because the 
motion cannot be considered as steady during the time interval 
2t:/lu. 




FIG. 5: Lyapunov exponents along the trajectory, for the solution 
drawn in Figure (T]-a). (a) negative real solution of equation ( I25t . 
(b)-(c) two other solutions, real parts in (b) and imaginary parts in 
(c). The real part of the two complex conjugate exponents crosses 
the zero real value at time tc — 2.3fo, here tc = 17.6475. 



where 



RESPONSE TO NOISE 



A{0o,vo) = 90 + [1 + e){l + liiivo)). 



(24) 



The eigenvalues of the matrix ( |23] ) are solutions of the equa- 
tion 



In view of forecasting catastrophes we have studied the re- 
sponse of the DRR system to an external white noise, as we 
did for the saddle-node bifurcation in the potential model 171]. 
We solved the set of stochastic equations 



u = v -1 + r]u.fu{t) 
v=^~^^(u+^{9 + \n{v)) 



= —V I 



+ Vvfv(t) 

+ {l + e)\n{v))+Tjgfe{t), 



(28) 



where /„, /„, fg are three independent noise functions with 
short memory time (of order St^qk) and rju,rj^,r]g are am- 
pHtudes smaller than 7^. We have calculated the correlation 
functions of a;(i), x standing for any function u, v, or 9, 

r^(i, t - t') =< x(t)x{t -t')> - < x{t) >< x{t -t') > , 

(29) 
which are functions depending on the delay t' but also on 
the time t, and the standard deviations 



a,it)=T,{t,t). 



(30) 



The results are very different from those of the potential 
model. First because here the standard deviation a^ of the re- 
sponse drastically increases with time before the event, much 
more than the velocity v{t) itself, that could be seen as a 
precursor. This contrasts with the potential case, where the 
standard deviation of the response grows nearly like v{t) and 
therefore cannot be used as a precursor The two functions 
cr^(i) and v{t), are drawn in Figure (|6}a) in In scale, showing 
that the growth of (t^ is three order of magnitude larger than 
the growth of v{t), on the time interval {tf — to, tf) , where 
the final time tf = tc — 0.1 corresponds to about two months 
before the event for typical earthquake case. 

Secondly the width of correlation function decreases very 
strongly during the large time interval ^ 2to before the event, 
see (|6}b). This could be understood as a sort of speed-up , con- 
trary to what happens in the potential model where the width 
of r„(i, t — t') becomes maximum shortly before tc, in agree- 
ment with the well-known critical slowing-down phenomena 
(or critical opalescence in spatial systems). 



correllation function 




maximum (critical slowing down) at time ^ tc — to (to be- 
ing the intermediate time scale in this system||7|]) then it drops 
abruptly before the catastrophe. For the DRR system we find 
that the correlation time of the response to noise (b-curve) de- 
creases very slowly until the singularity (speed-up). 

The slow decrease also lasts a time interval of order to, but 
here 



while it is 



FIG. 6: Response to noise (a)Comparaison of the increase of velocity 
and its standard deviation, in log scale,versus 'yt, until tc — 0.1, (b) 
correlation function Tu{t,t — t') as function of i' for three time i = 
tc — 2.3 (red), t = tc — 1 (orange), t — tc — 3O/7 (blue), for same 
parameters as in figure (TJ. 



to/Steqk ^ 7: 



to/Steqk ^ P^' 



(31) 



(32) 



for the saddle-node model. These two relations explain why 
the response to noise is changed much earlier in the DRR 
model than in the saddle-node model. Using parameters per- 
tinent for earthquake, the precursor time to is about few hours 
for the saddle-node model, while it is about few years for the 
DRR one. Both results are related to the linear stability of the 
SM since the memory time of the response to external noise 
is maximum when the temporal integral of the largest stabil- 
ity exponent is close to zero, that gives a "precursor time" of 
order to in both cases. 
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FIG. 7: Correlation time for (a) the potential model in ^, (b) the 
present model with 7 = 10^, .^ = 0.7,e = 0.72 



This striking result is summarized in Figure (|7|i which dis- 
plays the correlation time (mid-height width of correlation 
functions) of the fluctuating part of the displacement in both 
cases, along the trajectory. For the saddle-node model (a- 
curve) the correlation time increases slowly before it becomes 



POWER LAW FRICTION 

Finally it is worth questioning the relevance of this idea of 
transition for other systems. In other terms how "generic" is 
this kind of transition? One can say first that it persists (as a 



way of going from slow to fast and conversely for the full dy- 
namics) in the DRR equations by changing the numerical val- 
ues of the parameters. Another indication that it is a generic 
scenario is to check that it is still there in (slightly) modified 
equations. 
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FIG. 8: Solution of the DRR equations with ln(t;) friction (a), and 
their modified form with a power law friction (b) ( equations l |34t 
with a = 0.31) for the set of parameters 7 = 10, e = 1, ^ = 0.7. (c) 
\n{v) (red curve), and ga{v) (blue and green curves for a — 0.143 
and 0.333 respectively, (d) The ratio of the limit cycle periods (blue 
points) and amplitudes (green points) for the power law and ln(v) 
cases, as function of a 



The jacobian matrix is 







1 

-A{a) 






(36) 



where /^ = 



Aia) = 9o + il + e) 



6) ""<°+^)-\ and 



6'o(i),wo(i) are the dynamical variables along the trajectory, 
supposed to evolve very slowly with respect to the inverse of 
the imaginary part of any eigenvalue of the matrix(the Lya- 
punov exponents). We found that the eigenvalues of matrix 
(|36] | and (l23l l display similar behavior with one negative real 
exponent, and two others alternately real and complex. For 
example for 7 = 10"^, ^ = 0.7 and e = 0.78 we observe that 
the real part of the complex conjugate exponents cross zero at 
a time t^ — 2.45, which is very similar to the result obtained 
with the ln(w) friction term, cf figure ([5]). 

We have also studied the response to noise in the case of a 
power-law friction (equations ( |34l ) plus noise). The figure (|9]l 
displays a progressive speed-up, the correlation time decreas- 
ing during an interval of order few to before tc, as in the case 
of the ln(u) friction law (compare with Figure (|7]-b)). 

Correlation time 



We changed the logarithms dependence \n{v) in the DRR 
equations ^ and (|7]), by the function 



?a(w) 



(33) 



with a positive. Solving the modified system of equations 



FIG. 9: Half-height width of the correlation function, as a function 
oftimefora = 0.15,7 = 10^, e = 0.78,^ = 0.7 



ii ~ V — 1 

V ^ -j^ (u + ^{6 + g^iv)^ 

e = -V (d + {1 + e)g^iv)) , 



(34) 



we found stable periodic solutions, with slow and fast steps. 
More precisely this appears above a certain threshold (close 
to the Hopf bifurcation the limit cycle has ordinary behavior). 
The curves (a) and (b) of figure (HJ display slow-fast limit 
cycles obtained for the same parameter values, with the the 
ln(t;) friction and the power law friction, respectively (DRR 
equations and equations ( l34b respectively). In both cases we 
observe that the period of the limit cycle increases with | (for 
a given value of 7). Moreover in the power law case we find 
that the period of the limit cycle increases with a, see figure 
(c). This occurs of because the nonlinearity increases with 
the parameter a. The intermediate time (interval between the 
instant where v{t) = 1 and tc) is found to be nearly equal to 
unity. 



io-1, 



(35) 



SUMMARY 

We have shown on a model of stick-slip dynamics that RO 
can result from a breakdown of the adiabatic approximation. 
In such systems the transition from slow to fast motion is just 
the result of an unchecked acceleration of the slow motion, a 
rather natural effect that should appear beyond solid/solid phe- 
nomena. In the full DRR model, this acceleration is stopped 
by inertia, neglected in the adiabatic limit. In other contexts, 
the role of inertia could be played by feedback effects not 
taken into account in the equations of slow dynamics. Con- 
cerning the possibility of forecasting the fast event, and fol- 
lowing ideas presented in 17|], one finds that the response 
to an external source of noise is completely different in this 
model of stick-slip dynamics from what it is in the standard 
RO models. In the latter case the response to noise increases 
in amplitude before the transition and drifts to low frequen- 
cies (slowing-down). In stick-slip dynamics, the accelera- 
tion along the SM is also accompanied by a growth of the 



amplitude of the response, but it drifts to large frequencies 
(speeding-up). 
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